Loop quantum modified gravity and its cosmological application 
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A general nonperturvative loop quantization procedure for metric modified gravity is reviewed. As an exam- 
ple, this procedure is applied to scalar-tensor theories of gravity. The quantum kinematical framework of these 
theories is rigorously constructed. Both the Hamiltonian and master constraint operators are well defined and 
proposed to represent quantum dynamics of scalar-tensor theories. As an application to models, we set up the 
basic structure of loop quantum Brans-Dicke cosmology. The effective dynamical equations of loop quantum 
Brans-Dicke cosmology are also obtained, which lay a foundation for the phenomenological investigation to 
possible quantum gravity effects in cosmology. 
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I. INTRODUCTION 



As a background independent approach to quantize general relativity (GR), Loop quantum gravity(LQG) has been widely 
investigated in past 25 years CD-Hi- Recently, this non-perturbative loop quantization procedure has been generalized to the 
metric /(??) theories, JH 0] Brsns-Dicke theory ^ and scalar-tensor theories[8]. The fact that this background-independent 
, quantization method can be successfully extended to those modified theories of gravity relies on the key observation that these 
gravity theories can be reformulated into the connection dynamical formalism with a compact structure group. The purpose of 
this paper is to review how to get the connection dynamics of these modified gravity theories and how to quantize these theories 
by the nonperturbative loop quantization procedure. 

In fact, modified gravity theories have recently received increased attention in issues related to"dark matter", "dark energy" 
bi} and non-trivial tests on gravity beyond GR. Since 1998, a series of independent and accurate observations, including type la 
supernova, cosmic microwave background anisotropy, weak gravity lens, baryon oscillation, etc, implied that our universe now 
I , is undergoing a period of accelerated expansion! 9]. These results have caused the "dark energy" problem which is difficult to get 

■ a satisfactory interpretation within the framework of GR. Hence it is reasonable to consider the other possibility that GR is not a 
t^J- | valid theory of gravity on a galactic or cosmological scale. Besides the well-known /(ft) theories, a competing relativistic theory 

of gravity was proposed by Brans and Dicke in 1961 [ 10], which is apparently compatible with Mach's principle. To represent 
a varying "gravitational constant", a scalar field is non-minimally coupled to the metric in Brans-Dicke theory. To interpret the 
observational results, the Brans-Dicke theory was generalized by Bergmann [11] and Wagoner [12] to more general scalar-tensor 
7— I ' theories (STT). Moreover, scalar-tensor theories is also closely related to low energy effective actions of some string theory (see 
e.g. lfl3l4l5IO . Since it can naturally lead to cosmological acceleration in certain models (see e.g. lflq - [l9l0 . the scalar field in 
STT of gravity is expected to account for the "dark energy" problem. In addition, some models of STT of gravity may even 
account for the "dark matter" puzzle j20l - l22Tl . which was revealed by the observed rotation curve of galaxy clusters. 

■ There are infinite ways to modify GR. One may suspect which rules we should employed to do the modification. The decisive 
rule certainly comes from experiments. A large part of the non-trivial tests on gravity theory beyond GR is closely related to 
Einstein's equivalence principle (EEP) [23]. There are many local experiments existed in solar-system supporting EER which 
implies that gravity should be described by metric theories. Indeed, STT are a class of representative metric theories, which 
have been received most attention. That is why we us it as an example to demonstrate our general loop quantization procedure 
for metric modified gravity theories. It is also worth noting that both the metric f(fl) theories and Palatini fCR) theories are 
equivalent to certain special kinds of STT with the coupling parameter u>((/>) = and <jj(4>) — — \ respectively [24] Meanwhile the 
original Brans-Dicke theory is nothing but the particular case of constant co and vanishing potential of (p. There are also some 
other types of modified metric gravity theories proposed in recent years, such as Horava-Lifshitz theory 12511 and critical gravity 
rt26ll etc. Those theories are proposed based on the fact that GR is nonrenormalizable at perturbative level, while the introduction 
of higher order derivative terms might cure this problem. Thus it is quite interesting to see whether all those kind of metric 
theories of gravity could be quantized nonperturbatively. 

The following sections of this paper are organized as follows. We first introduce a general scheme of loop quantization for 
metric modified gravity in section [TT] Then we use scalar-tensor theories as an example to show how our general quantization 
procedure works. In section[lIl] we start with Hamiltonian analysis of STT. The coupling parameter of the STT naturally marks 
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off two sectors of the theory. Based on the resulted connection dynamical formalism of STT, we then quantize the STT by 
extending the nonperturbative quantization procedure of LQG in section [IV] in the way similar to loop quantum /(ft) gravity 
Nevertheless, the STT that we are considering are a much more general class of metric theories of gravity than metric 
fCR) theories. The Hamiltonian constraint operators in both sectors of the theory can be well defined. To avoid possible quantum 
anomaly, master constraint program of STT are given in[V] For cosmological application of above quantum gravity theories, we 
set up the basic structure of loop quantum Brans-Dicke cosmology and get its effective equations of motion in section[Vl] Finally 
some concluding remarks are given in the last section. Throughout the paper, we use Greek alphabet for spacetime indices, Latin 
alphabet a,b,c,..., for spatial indices, and i,j,k,..., for internal indices. 



H. GENERAL SCHEME 

In this section, we will first outline the general scheme of loop quantization for metric modified gravity l27ll . Especially, 
we are mainly focus on 4-dimensional metric theories of gravity which is consistent with Einstein's equivalent principle. The 
prerequisite is that the theory which is under consideration should have a well-defined geometrical dynamics, which means a 
Hamiltonian formalism with 3-metric h a \, as one of configuration variables, and in addition the constraint algebra of this theory 
is first-class (perhaps after solving some second-class constraints). Without loss of generality, we can assume that the classical 
phase space of this theory consists of conjugate pairs (h a i,, p ab ) and (<p B , n B ), where <p B could be a scalar, vector, tensor or spinor 
field. Then the quantization scheme has the following recipe. 

(i) To obtain a connection dynamical formalism, we first define a quantity K ab via 

K a b = [pdb ~ ^Phabj ■ (2.1) 

Then we enlarge the phase space by transforming to the triad formulation as 

(h ab ,p ab ) => (£? = VWUef, *i = K ab e h ). (2.2) 

Now we make a canonical transformation to connection formulation as: 

(E a j, H) => (E a j,Ai = r J a + jK ] a ), (2.3) 

and due to symmetric property of p a b we have K a [iE a ^ = 0. This will give us the Gaussian constraint, 2D a Ef = daE^+ei^A^E^ = 0. 
Then it is straightforward to write all the constraints in terms of the new variables, (ii) For loop quantization, we first represent 
the fields (<pB, n B ) via polymer-like representation, together with the LQG representation for the holonomy-flux algebra. Then 
the kinematical Hilbert space can be read as "74in := "HS, S'/t. All the basic operators and geometrical operators could be well 
defined in this Hilbert space. We can solve the Gaussian and diffeomorphism constraints as in standard LQG. Then we would 
get the gauge and diffeomorphism invariant Hilbert spaceas: 'K-^ n --» 'Kg "* Wolff- I n order to implement quantum dynamics, 
the Hamiltonian constraint operator may first be constructed at least in "Kg, although it usually could not be well defined in 
'Koiff- Then the master constraint operator can be constructed in 'Koiff by using the structure of the Hamiltonian operator, (iii) 
One may try to understand the physical Hilbert space by the direct integral decomposition of "Koiff with respect to the spectrum 
of the master constraint operator, (iv) One may also do certain semiclassical analysis in order to confirm the classical limits 
of the Hamiltonian and master constraint operators as well as the constraint algebra. The low energy physics is also expected 
in the analysis, (v) Finally, to complement above canonical approach, we can also try the covariant path integral (spinfoam) 
quantization. 

It should be noted that the last three steps are still open issues in the loop quantization of GR. Thus in the following sections, 
we will take scalar-tensor theories as an example to carry out the steps (i) and (ii) in the above scheme of loop quantization for 
metric modified gravity. 



III. HAMILTONIAN ANALYSIS OF SCALAR-TENSOR THEORIES 



The most general action of STT reads 



(3.1) 



where k - 8nG, % denotes the scalar curvature of spacetime metric g„ v , and the coupling parameter w(0) and potential £(<p) can 
be arbitrary functions of scalar field <p. Variations of the action ( 13.1b with respect to g ab and (f> respectively give us equations of 
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motion 



0G MV = V„V v - g^ncf, + ^-[(d^)d v <p - ^ V (V0) 2 ] - g^(<p), 



n + 



2a>(<f>) a>( 



-□<? 



(d^Wcf, + ^Mid^cf, - 2{'(<f>) = 0, 



7,2 ^ 



(3.2) 
(3.3) 



where a prime over a letter denotes a derivative with respect to the argument, V A , is the covariant derivative compatible with g MV 
and □ = g^VfjVy. By doing 3+1 decomposition of the spacetime, the four-dimensional (4d) scalar curvature can be expressed 
as 



% = K ah K ab -K 2 +R + -%=d M { ^n»K) - -?—d a ( yfhh ah d h N) 



(3.4) 



where K ab is the extrinsic curvature of a spatial hypersurface S, K = K a i,h ah , R denotes the scalar curvature of the 3-metric h a b 
induced on Z, rf is the unit normal of 2 and N is the lapse function. By Legendre transformation, the momenta conjugate to the 
dynamical variables h ab and <p are defined respectively as 



p ah = 



ft P \fh h ab 

i=- = -^WK ab - Kh ah ) - —(if, - N c d c <t>)], 

dh ab 2k N 
d£ ^h lv cj(<f>) 



7T = — ; = — — (k- 

d<p k N(f> 



• - N c d c 



where N c is the shift vector. The combination of the trace of Eq. (13. 5\ and Eq. d3.6t gives 

2kN 

(3 + 2w(0))(0 - N a d a <p) = —(4m - p). 



(3.5) 
(3.6) 

(3.7) 



It is easy to see from Eq. ( 13.7b that one extra constraint S — p — <pn — emerges when co((f>) = — I. Hence it is natural to mark 



off two sectors of the theories by co((f>) ¥= — § and w(</>) = — |. 



A. Sector of w(cf>) * -3/2 

In the case of (D(<p) + -3/2, the Hamiltonian of STT can be derived as a liner combination of constraints as 



HtoM = £d 3 x(N a V a +NH), 



where the smeared diffeomorphism and Hamiltonian constraints read respectively 
V(N) = £d 3 xN a V a = £d 3 xN a (-2D b (p ah ) + 7rd a (f>), 



H(N) 



1/ 



xNH 



d 3 xN 



■X 

By the symplectic structure 



2k 



PabP ah ~ \p 2 (p - <f>n) 



20(3 + 2cj) ) 



2k 



4> 



{h ab (x),p cd (y)} = ^6fd 3 (x,y), 
{<p{x)My)) = S 3 (x,y), 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



lengthy but straightforward calculations show that the constraints ( 13 .9b and ( D.lOl i comprise a first-class system similar to GR as: 



{V(N), V(N')} = V([N,N']), 
{H(M),V(N)} = -H(JL t M), 
{H(N),H(M)} = V(ND a M - MD a N). 



(3.12) 
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We can show that the above Hamiltonian formalism of STT is equivalent to their Lagrangian formalism|8]. 

To obtain the connection dynamical formalism of the STT, following the general scheme mentioned in last section, we define 



Lab uah 

K ab = 4>K ab + — - N c d c ab) = 4>K ab + -{<fm - p). 

2N (3 + 2w) Sh 



(3.13) 



Then we can introduce new canonical pairs (E" = -\fhe",K' a = K a i,e h ), where e" is the triad such that h a be1e b = tfy. Now the 
symplectic structure (13. lit can be got from the following Poisson brackets: 



{E a j{x),E b k {y)} = {K J a {x),K k h {y)) = 0, 
{kl{x),E b k {y)} = Kd b a 6 i k 6{x,y). 

Note that since we have K ah = K ba , there have an additional constraint: 

G jk = K a[j E a k] = 0. 

Now we can make a second canonical transformation via defining: 

K=K + yR i a , 



(3.14) 



(3.15) 



(3.16) 



where T' a is the spin connection determined by E", and y is a nonzero real number. It is easy to see that our new variable 
A J a coincides with the Ashtekar-Barbero connection l28l l29ll in the special case = 1. The Poisson brackets among the new 
variables read: 



{A J a {x),E b {y)} = yK6 b 6 J k 6{x,y), 
{A^xlAiiy)) = 0, [E a j{x),E b k (y)} = 0. 



(3.17) 



Now the phase space of the STT consists of conjugate pairs {A' a ,E b ) and (<p,n). Combining Eq. (l3. 15b with the compatibility 
condition: 



d a E f + e ijk TiE ak = 



the standard Gaussian constraint can be obtained as 



j Tpak 



Qi = @ a E? = d a E? + e ijk A J a E 



(3.18) 



(3.19) 



which justifies A' a as an sM(2)-connection. It is worth noting that, had we let y = ±i, the (anti-)self-dual complex connection 
formalism would be obtained. The original vector and Hamiltonian constraints can be written in terms of new variables up to 
Gaussian constraint respectively as 



1 

icy 



V a = —F^+nda 



(3.20) 



H = f K 



Kb-(y 2 + ^ s ^ R 'a n Kb 



SjuE a k E b 
a/A 



'{K^f | 2 {K' a E«)n | K 



(3 + 2u>(qb)) 



K 2 qb\Th 



1 

+ - 

K 



(4<f>) 



2(f> 



yfh(D a (p)D a (p + yfhD a D a <p + Vft£(0) 



where F' ab = 2d[ a A' b , + e' kl A a A h is the curvature of A' a . The total Hamiltonian can be expressed as a linear combination 



H tota i = jA i @i+N a V a + NH. 



(3.21) 



(3.22) 



It is easy to check that the smeared Gaussian constraint, £?(A) := ^d i xA'{x)@i{x), generates SU{2) gauge transformations on 
the phase space, while the smeared constraint, 'V(iV) := L d 3 xN a {V a - A' a Qi), generates spatial diffeomorphism transformations 
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on the phase space. Together with the smeared Hamiltonian constraint H(N) - L d 3 xNH, we can show that the constraints 
algebra has the following form: 

{0(A), 0(A')} = *0([A,A']), (3.23) 

[0(A), <V(lk} = -eU-jA,), (3.24) 

[0(A), H(N)} = 0, (3.25) 

{'V(N),^V(N')} = <V([N,~N']), (3.26) 

{fV(N),H(M)} = H(£^M), (3.27) 
{H(N),H(M)} = < V(ND"M - MD a N) 

+ g{(Nd a M-Md a N)h ab A b ) 
[E a D a N, E b D b MY 



Kh 



Si 



2 [E a D a ($N),E»D b (4>M)}' r 
- r ; m- (3.28) 

Here Eqs. (13.2313. 27b can be understand by the geometrical interpretations of Q{A) and 'V(N). The detail calculation on the Pois- 
son bracket ( 13.28b between the two smeared Hamiltonian constraints can be seen in the Appendix of 08J . Hence the constraints 
are of first class. Furthermore, the constraint algebra of GR can be recovered for the special case when </> — 1. To summarize, the 
STT of gravity in the sector cl>(<P) + -3/2 have already been cast into the si<(2)-connection dynamical formalism. The resulted 
Hamiltonian structure is quite similar to metric f(R) theories 16J. 

B. Sector of cj(cf>) = -3/2 

In the sector of a>((p) = -3/2, Eq. (13.71 i implies that there is an extra primary constraint S = 0, which we call "conformal" 
constraint. Hence, as pointed out in Ref. lBOIl . the total Hamiltonian in this case can be expressed as a liner combination of 
constraints as 



Hiatal - J ' 



d 3 x(N a V a +NH + AS ), (3.29) 

where the smeared diffeomorphism constraint V(N) is as same as ( 13.91 l. while the Hamiltonian and conformal constraints read 
respectively 



H(N) = I d 3 xNH 



I 



2k 


(Pabp ab - 


{P2 \ 


yfh 


{ 


) 


-\ 


d 3 xA(p - 


<fm). 



^fh 3 
+ —(-$R - —{D a <t>)D a $ + 2D a D a 4> + 2f(«) 
2k 2<p 



(3.30) 
(3.31) 



With the help of symplectic structure ( 13.111 ). straightforward calculations show that 



{H(M),V(N)} = -H{JL t M), {S(A),V(N)} = -S(X 1 z<l), (3.32) 
{H(N),H(M)} = V(ND a M - MD a N) + S(^-(ND a M - MD a N)), (3.33) 



IS (A). 



AM C r 

H{M)} = H{—) + J NAyfh(-2t(<l>) + </>?(</>)). (3.34) 



The Poisson bracket (13.34b implies that, we have to impose a secondary constraint in order to maintain the constraints S and H 
in the time evolution as 

-2#0) + ^'(0) = O. (3.35) 



It is clear that this constraint is second-class, and hence one has to solve it. Here we consider the vacuum case where the solutions 
of Eq. ( 13.35b are either £(</>) = or£(0) = C<p 2 , where C is certain dimensional constant. Thus the consistency condition strongly 
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restricted the feasible STT in this sector to only two theories. As pointed out in Ref.[24], for these two theories, the action ( 13.11 ) 
is invariant under the following conformal transformation: 

g,j V -> e A g IJ y, 4> e~ A 4>. (3.36) 

Thus, besides diffeomorphism invariance, those two theories are also conformally invariant. Now in the resulted Hamiltonian 
formalism the constraints {V, H,S) form a first-class system. The following transformations on the phase space are generated by 
the conformal constraint 

{h ab ,S(A)} = Ah ab , {P ab ,S(A)} = -AP ab , (3.37) 
{cf>,S(A)} = -A<f>, {n,S{A)) = An. (3.38) 

It is clear that the above transformations coincide with those of spacetime conformal transformations (13.36b . Hence all constraints 
have clear physical meaning. Now the physical degrees of freedom of this special kind of STT are equal to those of GR because 
of the extra conformal constraint (13.3 It . 

The connection-dynamical formalism for STT in this sector can also be obtained by the canonical transformations to the new 
variables (13.16b . Then the total Hamiltonian can be expressed again as a linear combination 

a tota i = f KQi + N a V„ + NH + AS, (3.39) 



where the Gaussian and diffeomorphism constraints keep the same form as Eqs. ( 13.191 1 and (13.20b , while the Hamiltonian and 
the conformal constraints read respectively 



e Jkl E l E l 



K 

Now the constraints algebra in the connection formalism is closed as 



+ i[ - y/h(D a <l>)D a (f> + VhD a D a (f> + V^(0)], (3.40) 
S = -K' a E" - n<p. (3.41) 



{0(A), H(N)} = 0, (3.42) 
{0(A), S (A)} = 0, (3.43) 

{S(A),H(M)} = fl(^), (3.44) 



{H(N),H(M)} = [ r V(ND a M - MD a N) 



+ S(—(ND a M - MD a N)) 

+ Q({Nd a M-Md a N)h ab A,) 
[E a D a N,E b D b M] i 

, [E a D a (d>N), E b D h {<bM)Y , 
- y- — ^—^e^. (3.45) 

It is easy to see that the Poisson brackets among the other constraints weakly equal to zero. Hence we have cast STT of gravity 
with iL>((p) + -3/2 into the s«(2)-connection dynamical formalism. 



IV. LOOP QUANTIZATION OF SCALAR-TENSOR THEORIES 



Based on the connection dynamical formalism obtained in last section, the nonperturbative loop quantization procedure can 
be naturally extended to the STT. The kinematical structure of STT is as same as that of LQG coupled with a scalar field and 
fCR) theories |5j,|6[]. The kinematical Hilbert space of the system is a direct product of the Hilbert space of geometry part and 
that of scalar field, 'H^n := HE^ ® *7fff n . An orthonormal basis of this Hilbert space is the so called spin-scalar-network basis, 
T a x(A, 0) = T a (A) ® Tx(<p), over some graph aUlcE, where a and X consist of finite number of curves and points in 2 
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respectively. The basic operators of quantum STT are the quantum analogue of holonomies h e (A) — !Pexp J e A a of a connection 
along edges e c S, densitized triads E(S,f) := J s e a b c E"f smeared over 2-surfaces, point holonomies Ua = exp(//t</>(x)) ll3~lll . 
and scalar momenta n(R) := J R d 3 x7r(x) smeared on 3-dimensional regions. It is worth noting that the spatial geometric operator, 
such as the area[32] , the volume[33] and the lengthily l35ll operators, are still valid in of quantum STT. As in LQG, 
it is natural to promote the Gaussian constraint ^(A) as a well-defined operator^ 0]. Then it's kernel is the internal gauge 
invariant Hilbert space 'Hq with gauge invariant spin-scalar-network basis. Since the diffeomorphisms of X act covariantly on 
the cylindrical functions in "Hq, the group averaging technique can be employed to solve the diffeomorphism constraint yl|4l|6|]. 
Hence the desired diffeomorphism and gauge invariant Hilbert space < Hmff for the STT can also be obtained|0,[^]. 



A. Sector of cj(<f>) * -3/2 

While the kinematical framework of LQG and polymer-like scalar field have been straight-forwardly extended to the STT, 
the nontrivial task in the case of u>((p) + -3/2 is to implement the Hamiltonian constraint ( 13.211 ) at quantum level. In order 
to compare the Hamiltonian constraint of STT with that of metric fCR) theories in connection formalism, we write Eq. ( 13.211 ) 
as H(N) = 2f=i Hj in regular order. It is easy to see that the terms Hi,H2,H-i,H% just keep the same form as those in f{K) 
theories (see Eq.(39) in Ref.|01), and the H^,Hn,Hi terms are also similar to the corresponding terms in fCR) theories. Here 
the differences are only reflected by the coefficients as certain functions of </>. Now we come to the completely new term, 
H§ = j^d^xN^^- Vh(D a (p)D a <p. This term is somehow like the kinetic term of a Klein-Gordon field which was dealt with in 

Ref.fUl]. We can introduce the well-defined operators (f>, <p~ [ as in Ref. jj^l . It is reasonable to believe that function a>(<p) can 
also be quantized [6]. By the same regularization techniques as in Refs.|g|3g|, we triangulate X in adaptation to some graph a 
underling a cylindrical function in fi^in and reexpress connections by holonomies. The corresponding regulated operator acts on 
a basis vector T a jc over some graph a U X as 

H b T a , x = lim (f> (v)# e (v-v) 

e^O 3 y \lA()f ylTl) K 



EM 

vea(y) v ' v(A)=v 
X [U Ao {(p(t SL (A v ))) ~ U Ao {(p(s SL (A v) ))] 

x Tr(T4 M(Av) [^;2 (Av) , (VuO V4 ]K(aAKU,)> (Vu;.) V4 ]) 

v'ea(v) ( ' v(A')=v' 

x TrCT^,,^)^; 1 ^,), {Vu^ VA ]K(AAKl(A,y (%) 3/4 D ■ T a ,x- (4.1) 

We refer to |6] for the meaning of notations in Eq.( |4.1t . It is easy to see that the action of on T ar x is graph changing. 
It adds a finite number of vertices at f(i/(v)) = e for edges ej(t) starting from each high-valent vertex of a. As a result, 
the family of operators H^(N) fails to be weakly convergent when s — > 0. However, due to the diffeomorphism covariant 
properties of the triangulation, the limit operator can be well defined via the so-called uniform Rovelli-Smolin topology induced 
by diffeomorphism-invariant states $>Diff as: 

® Dif f(H 6 ■ T a , x ) = lim((D D , // |// 6 E |r ff , x ). (4.2) 
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It is obviously that the limit is independent of e. Hence both the regulators e and s can be removed. We then have 

KA)=v(A')=v 
X [^ o (0(f. Si (A))) - U Ao ((p(s sdA) ))] 
x Tr(r4 M (A)[/^ (A) , (^) 3/4 ]^(A)[^ (A) , (K) 3/4 ]) 
x e( S/ i / ^)e / ^^- 1 (0(^ ;(AO )) 

X [^ o (0(f S;( A'))) - C^M-Va')))] 

X Tr^XAO^W (Vv) 3/4 ]^ SJC (A0^;](A0' (^) 3/4 D ' ^ ( 4 -3) 

In order to simplify the expression, we introduce f((f>) = 3+2 ^ w f° r me other terms containing it in H(N), which can also be 
promoted to a well-defined operator f ((/>). Hence, the terms H3, H4 and H5 can be quantized as 

A T v 4 W0/WV)) 3L-1, n 

x [H E (l),ft v ][HE(l),Jp v ]-T a , x , 

(4.4) 



x J] Tr^A)^,,!]) 

KA)=v(A')=v 

LAW 



x Tr(T ( ft SM(A) [ft;2 ( A). (n) 3/4 ]/j. !N (A) [h;l (A) , (v v ) V4 ]) 

x ^SjSk)^ 

x Trih^ih;^, {%f' 2 ]h sjm [h- s ] m , (V v ) 1 ' 2 ] 

x ^(Aot^AO'^v) 172 ])-^, (4.5) 



h t V 2 1 VjV(v)/(0(v)) ^ 

X J] ^SjSkW* 

v(A)=v(A')=v 

x Tr(^. !/(A) [^; (A) , (y v ) 1/2 ]/*. !j( A)[/*;; (A) , (v,) 1/2 i 

x ^(A)[^ (A) ,(K) 1/2 ]) 

x e{s L s M s N )^ MN 

x TrC^Aot^V). (K) 1/2 ]^„ ( ao[^ ( ao> (^') 1/2 ] 

x Km\Kl(KY^n-T a , x . (4.6) 
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While the H-j and //g terms keep the same form as in fCR) theory, which read respectively 

= y 2 7 /V(v) 

7 "' X 3r 2 ^o(^) 2 ^(v) 

x Z x ' Z ^/*)« ,/jr 

e(0)=v v(A)=v 
X f>^(^ S/ (A)))[^o(0fe/(A))) ~ #<lo(0O*,(A)))] 

X T*t,S,^)[^ A) ,(P' v ) 1/2 ] 

x ^MA)[^ (A )'(K) 1/2 ])-^,x, (4.7) 



Hs ■ = - V 7Y(v)|(0(v))f v ■ r«, z . (4.8) 

veV(a) 

Here the action of the volume operator V on a spin-network basis vector T a (A) over a graph a can be factorized as 

V-T a = ^ t/,,-r„. (4.9) 

Therefore, the total Hamiltonian constraint in this sector has been quantized as a well-defined operator H(N) = 2)f=i @i m ^kin- 
It is easy to see that H(N) is internal gauge invariant and diffeomorphism covariant. Hence it is at least well defined in the 
gauge invariant Hilbert space Iri^. However, it is difficult to define H(N) directly on Hoiff- Moreover, as in f(R) theories, the 
constraint algebra of STT do not form a Lie algebra. This might lead to quantum anomaly after quantization. 



B. Sector of aj(cf>) = -3/2 



In the special case of w(</>) = -3/2, the smeared version S (A) of the extra conformal constraint ( 13.41b has to be promoted as a 
well-defined operator. Note that both <p and n(R) are already well-defined operators. We can use the following classical identity 
to quantize the conformal constraint S (A), 



K 



= Jd*xkiEf = y-HH E (l),V} 



(4.10) 



where the Euclidean scalar constraint H E (l) by definition was: 



H E {\) 



2k L 



d\F j 



apb 



(4.11) 



Both H E and the volume V have been quantized in LQG. Now it is easy to promote S (A) as a well-defined operator, and its 
action on a given basis vector T a ^ e 'Hidn reads 



n.X 



KveV(a) 



(4.12) 



It is clear that S (A) is internal gauge invariant, diffeomorphism covariant and graph-changing. Hence it is well defined in < Hq. 
The Hamiltonian constraint operator in this sector is similar to that in the sector of w(</>) + -3/2. The difference is that now 
= -3/2. Hence we write Eq. d340l as H(N) = £f =1 Hi in regular order. It is easy to see that the terms H\,H2, H4, H5 just 
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keep the same form as those in last subsection, while the quantized version of H3 is 

A t - V 2 1( W(v) % _ lf 

x J] e(5 i 5 M ^)e iMJV ^ 1 (0(5 Si (A))) 

KA)=v(A')=v 
X [U Ao (<f>(ts L (A))) - Ua (4>(Ssl(A)))] 

x Tr(r4 M (A)[/^ (A) , (<^v) 3/4 ]^[^, (K) 3/4 ]) 

X [^ O (0(?„(A'))) - ^(0(*fc(A')))l 

X TltT4 j( A<))[^(A<)' (Vv) 3/4 ]^ Sr CA0^^(A0' (^) 3/4 D ' ^ ( 4 - 13 ) 

Hence the total Hamiltonian constraint operator H(N) = 2/=i "i now is also we U defined in Hq. 



V. MASTER CONSTRAINT 



In order to find the physical Hilbert space and avoid possible quantum anomaly, master constraint programme was first intro- 
duced into LQG by Thiemann in his seminal paper l37ll . The master constraint can be employed to implement the Hamiltonian 
constraint. This programme can also be generalized to the above quantum STT. 



A. Sector of to(<f>) j= -3/2 
In the sector of a>((p) + -3/2, by definition, the master constraint of the STT classically reads 

KA 1 Cj3 I^WI 2 ,c n 

M:= 2J/ X ^' (5A) 

where the expression of Hamiltonian constraint H(x) is given by Eq. (13. 2U . The master constraint can be regulated by a 
point-splitting strategy ll38ll as: 

M= l -f£y[d^( X -y)^=^=. (5.2) 



Introducing a partition P of the 3-manifold E into cells C, we can get an operator Hg„ acting on the internal gauge-invariant 
spin-scalar-network basis T s c in "Kg via a state-dependent triangulation, 

H s c, a ■ T s , c = J] Xc{v)H* ■ T s , c , (5.3) 
where xc is the characteristic function over C, a denotes the underlying graph of the spin-network state T s , and 

#v= z #st + Z% (5 - 4) 

KA)=v i=3 

with 

ftE 16/Wv)) ,. 

x [^(1),(^) 1/4 ][#*(1),(Pb.) 1/4 ], (5.5) 
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y 2 18 /(^(V)) , ) 

x e(s L s M s N )e LMN Tr(Trh SM(A) [h^ (A) , (V m ) 1/2 ] 

x W^WW 1 ' 2 ]) 

x e{s I sjs K )e lJK TxOi Sim [h-] m , (V ut ) 112 ] 

x W)[^),(^) 1/2 ] 

x Km&'s^Ay&ui) 1 ' 2 ]), (5.6) 



v 2 20 Kf(cf>(v)) 
KA)=v(AO=v(X)=v 3 V(^) 6 ^ 2 (V) 

x e(^ y ^) e "*Tr(^ ;(A) [^/ (A) , (V^ ) 1/4 ] 

x W^W^) 1/2 ]) 

x e( iL . M ^)6 iMA 'Tr(^. Si(A0 [^ i 1 (A , ) , (Vf/,0 174 ] 

x V(A')[^ ( ao'(^) 1/2 ] 

x W)fc w W(^) 1/2 ])' ( 5J ) 



#6,v 



y 2' 5 &(0) . 



v(A)=v(A')=v 

x e( iL . SM ^)e iMA 'j7; o 1 (0(^ i(A) )) 

x [^ o (0(f Si( A))) - U Ao {(p(s SL(A) ))] 

x TrCr^^^t^, (W' 2 ]^^, (K) 3/4 ]) 

X [^(^(^(AO)) - #4,(0C*«(AO))] 

x Tr(T^ [^ , (^v) 1/2 ]W)[^). (^) 3/4 ])' ( 5 - 8 > 



7 ''' 3r 2 //io(^) 2 ^£(v) 
x E*Z 

e(0)=v v(A)=v 
X [t>^(0(^ ;( A))) - ^(0(j ft (A)))] 

x Tr(r,/i. !j(A) [/i;; (A) ,(y (/ , e ) I/4 ] 

x ^(A)[^ (A) ,(^) 1/4 ]), (5.9) 

% = i#(0(v))^"; (5.10) 

Here ^^r v keeps the same form as the corresponding terms in @]. Note that the family of operators H E Ca are cylindrically 

consistent up to diffeomorphism. Thus the inductive limit operator He is densely defined in < Hc by the uniform Rovelli-Smolin 
topology. Hence we could define master constraint operator M. acting on a diffeomorphism invariant state as 

<A<bmff)T,< : = lim <D w// [iV H S C (H^T S , C ]. (5.11) 
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Note that although the quantitative actions are different, our construction of M is qualitatively similar to those in JHH^)]. Similar 
methods in (a, |36J] can be used here to prove that M. is a positive and symmetric operator in I^Diff- Hence it admits a unique 
self-adjoint Friedrichs extension. It is then possible obtaining the physical Hilbert space of the quantum STT of this sector by 
the direct integral decomposition of 'Huij j' with respect to M. 



B. Sector of = -3/2 



In the case of a>(<p) = -3/2, both the implementation of the Hamiltonian constraint and the implementation of the conformal 
constraint need to employ the master constraint programme. We then define the master constraint for this sector as 

M:= i r dix \Hi X t + \s {x f 

where the expressions of Hamiltonian constraint H(x) and the conformal constraint Six) are given by Eqs. (13.401 1 and (13.411 ) 
respectively. It is clear that 

W( = o //(/V) = and S(A) = VN(x),A(x). (5.13) 
Now the constraints form a Lie algebra. The master constraint can be regulated by a point-splitting strategy as: 

M = o d y d x Xe(x-y) . . (5.14) 



Introducing a partition V of the 3-manifold 2 into cells C, we get an operator acting on spin-scalar-network basis T sx in 
<Hc by a state-dependent triangulation as Eq. ( 15.31 ). Here, note that Hf, has less terms than in Eq. (15.31 ) as 

a e = 2 #st+Z% (5 - 15) 



v(A)=v i=3 



where 

- _ y 2 14 ,_j 

X eiSLSMSN)^ 0^(<f>(S Si (A))) 

x [^(^(^(A))) - U Ao ((p(s sdA) ))] 

x Trfo^t^. (K) 1/2 ]^ !W( A)[^; (A) , (K) 3/4 ]) 

x eisjSjSK^U^icpiss,^))) 

x [UaM 1 ^')))- UaM s s,(A')))] 

X TT(Tih Sjm [h-j m , (K) 1/2 ]^,(A')[^(AT (K) 3/4 ]), (5.16) 

and H^ v and Hf keep the same form as the corresponding terms in the sector of a>(4>) + -3/2. The operator corresponding to 
the conformal constraint can be defined in a similar way, 

S E c, a -T s , c = J] Xc(v)S E v -T s , c , (5.17) 

veV(a) 

where 

St = §i* +S iv> (5 - 18) 

with 

S lv = 3/2 Lj ^ £ (D,(^) 1/2 ], (5-19) 
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3r 3 (^) 3 £(v) 
x eCi/i/iri^'TrC^E^, (V^) 1/2 ] 
x ^,(A)[^; (A) ,(V^) 1/2 ] 

x V(A)[^ ( a).(^) 1/2 ])- (5-20) 

Note that the family of operators H e Ca and SJ a are cylindrically consistent up to diffeomorphism. Hence the inductive limit 

operator He and S c can be densely defined in < Hg by the uniform Rovelli- Smolin topology. Thus we could define master 
constraint operator M. acting on diffeomorphism invariant states as 

(M$Dtff)T s ,c = Um Y. {h s c (H^ + S C (S^) T s , c ]. (5.21) 

Similarly, we can prove that M. is a positive and symmetric operator in 'HdiJ] and hence admits a unique self-adjoint Friedrichs 
extension|0[3^]. Hence it is also possible to obtain the physical Hilbert space of the quantum STT in this special case by the 
direct integral decomposition of 'Hoiff with respect to the spectrum of AT 



VI. COSMOLOGICAL APPLICATION OF QUANTUM BRANS-DIKCE THEORY 



For cosmological application of above loop quantum STT, in this section, we will set up the basic structure of loop quantum 
Brans-Dicke cosmology and get its effective equations of motion II39I1 . For simplicity, we only consider the spatially flat (k = 0) 
homogeneous and isotropic universe in Brans-Dicke theory. Recall that the original Brans-Dicke theory is the particular case of 
STT with constant u> and vanishing potential of <p. Thus the Hamiltonian constraint of Brans-Dicke theory reads 



H = — 

2k 



b 

k^l 



3 + 2u 



+ ■^-^Jh(D a <p)D a <p+-TjhD„D a q 

lK(f> K 



(6.1) 



In the cosmological model, classically the metric of spacetime can be written as the following Friedman-Robertson- Walker 
(FRW) formalism, 



ds 2 = -dt 2 + a 2 (t)[dr 2 + r 2 (d0 2 + sin 2 0dcf> 2 )] 
where a is the scale factor. The classical Friedman equation of Brans-Dicke cosmology reads 



a cp 
a 2<p 



2w + 3 
12 



- + 



(6.2) 



(6.3) 



Our task is to quantize this model by the loop quantization method and find the quantum dynamical equation as well as its 
effective expression. 



A. Loop quantum Brans-Dicke cosmology 

Since the space of our cosmological model is infinite, we introduce an "elemental cell" *V and restrict all integral to *V. The 
homogeneity of the universe guarantee that the whole space information is reflected in this elemental cell. Now we choose a 
fiducial Euclidean metric "q a b and introduce a pair of fiducial orthnormal triad and co-triad as (°ef, °oj' a ) respectively such that 
"lab = <^>'a 0<x >h- F° r simplicity, we let the elemental cell *V be a cubic measured by our fiducial metric and denotes its volume as 
V . Because our FRW metric is spatially flat, we have T' a = and hence A' a = yK' a . Via fixing the degrees of freedom of local 
gauge and diffeomorphism, we finally obtain the connection and densitized triad by symmetrical reduction as l40ll : 



N a = cV^X, E) = pvj ^S)°e), 



(6.4) 
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where c, p are only functions of t. Hence the phase space of the cosmological model consists of conjugate pairs (c, p) and (<p, n). 
The basic Poisson brackets between them can be simply read as 

{c,p} = -y, 

{•M) = 1. (6.5) 

Note that by the symmetric reduction, the Gaussian and diffeomorphism constraints are satisfied automatically. Also in our 
homogeneous model, the last two spatial derivative terms in the Hamiltonian constraint ( 16.11 ) can be neglected. Hence we only 
need to consider the first five terms in ( 16. 11 1. The reduced Hamiltonian in the cosmological model reads 

3c 2 VLpI k ,3cp 2 
H = + -(— +n<f>y. (6.6) 

y ^ (3 + 2w)#7p K y 

To quantize the cosmological model, we first need to construct the quantum kinematic of Brans-Dicke cosmology by mim- 
icking the loop quantum STT This is the so-called polymer-like quantization. The kinematic Hilbert space for the geometry 
part can be defined as < K,^ 1 „ := L 2 (RBohr, dpH), where Rsohr and dpH are respectively the Bohr compactification of the real line 
and Haar measure on it Boll . On the other hand, for convenience we choose Schrodinger representation for the scalar field |4l||. 
Thus the kinematic Hilbert space for the scalar field part is defined as in usual quantum mechanics, < K^ n := L 2 (R,dp). Hence 
the whole Hilbert space is a direct product, < Hl°™ 1 - ( H^ B ® "H^L- Now let \fx) be the eigenstates of p in the kinematic Hilbert 
space ftf? such that 

m = *^,\p). (6.7) 
6 

It turns out that those states satisfy the following orthonormal condition 

(w\Hj) = , (6.8) 

where 6 MiJJ . is the Kronecker delta function rather than the Dirac distribution. For the convenience of studying quantum dynamics, 
we define new variables 



v :- 



2 V3sgn(p)/T 3 , b := pc, (6.9) 



where sgn(p) is the sign function for p and p = with A = 4 yf^nyf 2 being a minimum nonzero eigenvalue of the area 
operator J42j]. They also form a pair of conjugate variables as 



{b,v} = \- (6.10) 
n 

It turns out that the eigenstates of v also contribute an orthonormal basis in < H^.. We denote \</>, v) as the orthogonal basis for the 
whole Hilbert space H'^'. 

Now we come to the quantum dynamics. We treat the first two terms of Hamiltonian constraint (16.1b in the same way as in 
standard LQC l43ll . Hence, the first two terms of the Hamiltonian constraint act on a quantum state *P(v, <p) e 'K^ 1 as 

(#! + H 2 W(v, <t>) = \ CAOWv + 4, 0) + / (v)T(v, 0) + /_(v)T(v - 4, <j>)) , (6.11) 



where 



V3A | 

/-(v)=/ + (v-4), /o(v) = -/ + (v)-/_(v). (6.12) 



/ + (v)=T7T2 Iv + 3|-|v+l||v + 2|, 



Then we turn to the third term = =-4 — ° V . Due to spatial flatness, we have K'E a = -A'E a . In the cosmological model, 
this term can be reduced by 



1 . „ 3 3Khbv 

—A 'E w — cp = 

y a ' y 1 4 



~A' a E? w -cp = ——. (6.13) 
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Because we use polymer representation for geometry, their is no quantum operator corresponding to connection c as in standard 
LQC [4lJ . Hence we have to replace the connection by holonomy to get a well-defined operator. It turns out that the term 
can be quantized, and its action on a quantum state reads l39ll 

2 a/3" / 3 fi \^ 
A 3 *(0,v) = — — -r sm(b)\v\ S mQ>m<t>,v) 
MA)§ \ 4 / 

\Fx I \^ 

1K T [|v + 2|^(^, v + 4) - 2^1^(0, v) + |v - 21^(0, v - 4)] , (6.14) 



where we set j3 = 3 + 2<w. Similarly, the fourth term H4 = —j^fc- can also be quantized, and its action on a wave function 
reads 

2V3W37A 

H 4 W((p,v) = r — 2sgn(p) sin(fc)^(0, v) 

j8(A)i \ 4 / 

= r — Ssgn(p)[ — — ]. (6.15) 

J3(A)i \ 4 / 50 30 

The last term H5 = ^= can be quantized as 

A A) I 



where 



2 V3/c , d 2 ^, v) 
-^W 2 g(v)0 ; , (6.16) 

/3(A) 2 30 2 



fi(v) = (?) 3 |v|||v+l| 1/3 -|v-ir /3 | 3 . (6.17) 



"2 

The total Hamiltonian constraint equation of loop quantum Brans-Dicke cosmology reads 



(J]8 i mcf>,v) = 0. (6.18) 



B. Effective equation 

To study the effective theory of loop quantum Brans-Dicke cosmology, we also want to know the effect of matter fields on the 
dynamical evolution. Hence we include an extra massless scalar matter field tp into Brans-Dicke cosmology. Then classically 
the total Hamiltonian constraint of the model reads 

3c 2 VH k ,3cp - V\ 

H= 7— + r( — + 7r 0) + — r> ( 6 - 19 ) 

y- K <f> (3 + 2co)cf>\p\-- K r 2\ P \-2 

where p v is the momentum conjugate to <p. The effective description of LQC is a delicate and topical issue since it may relate 
the quantum gravity effects to low-energy physics. The effective equations of LQC are being studied from both canonical 
perspective [44-47] and path integral perspective[48-52]. Since the key element in the polymer-like quantization of previous 
subsection is to take holonomies rather than connections as basic variables, a heuristic and simple way to get the effective 
equations is to do the replacement c — > sm( ^ c) or b — > sin b. Under this replacement, the effective version of Hamiltonian 
constraint d6.19t takes the form 

3 sin 2 (fie) J\p\ k 3sin(Z/c)n , 3 

H = T,-Y + r( - + ^) + Wp, (6-20) 



16 



where p = by definition is the matter density. It is worth noting that the effective Hamiltonian ( 16.201 can also be derived by 

z\p\ 



a path integral formalism [39]. Then the canonical equations of motion read 

^3 sin(fic)p 
pKy 



p = — sm(pc) cos(pc) r sgn(p) 



7# 

2k 3 sm(flc)p 
-( — +nd 



P4>W 



+ 7T(p) COS(jUC 



flp\* ^7 

In the above calculation, the Poisson brackets (16.51 ) were used. The Combination of equations d6.2U and ( 16.22b gives 



P_ _0_ 

2p 20 



sgn(p) 

y<f>p Vipi 

sgn(p) 



sin(/2c) cos(pc) 



k 3sin(pc)p 

-( — z + - co&ific)) 



P«y 



sin(pc) cos(pc) + —(1 - cos(pc)) 
y0 VA 20 

On the other hand, from effective Hamiltonian constraint ( 16.201 we can get 

3 sin V) ffi 

+ hp = 0, 



which implies 



Kj 2 (pA 4/ap 



si nV) = 



(6.21) 
(6.22) 



(6.23) 



(6.24) 



(6.25) 



where p c = 



V3 



y 2 A/f 327r 2 G 2 y- 1 ?! 



3^ and p e yy = ^ + 0p. Taking account of Eq. d6.25l >, we can rewrite Eq. d6.23t as 

2 



a <p 
a 20 



1 \hPeffi l ) + T7(! - a/ 1 ) 

V 3 p c 20 V Pc 



(6.26) 



This is the effective Friedmann equation of Brans-Dicke cosmology, which contains important quantum correction terms. In 
addition, we can show that for a contracting universe, p e ff monotonically increase while v decreases(39J]. Thus it is easy to see 
from Eq. (16.26b that, when p e jf approaches p c , one gets cos(/ic) = 1 - ^ = 0. Then from Eq. d6.2U , we can obtain p = 0. This 
implies a quantum bounce would happen at that point for a contracting universe. 

We end up this section with several remarks. First, when 0=1, because of p e ff = p and = we would return to the 
well-known effective Friedmann equation of LQC ||43j |45[] as 



ia\ 2 k p 

- = TPd--). 

V a I 3 p c 

Second, when p e ff «: p c , we can omit terms in Eq. ( 16.26b to get the classical limit of this equation as 



(6.27) 



a \ 1 k k B(f> 2 

= ^ i Kp 
120 2 30' 



(6.28) 



which is nothing but the classical Friedmann equation of Brans-Dicke cosmology. Hence the effective theory has correct classical 
limit. 



VII. CONCLUSION AND OUTLOOK 



Modified gravity has received increased attention in issues of "dark matter", "dark energy" and nontrivial tests on gravity 
beyond GR. Some kinds of modified gravity theories have also become popular in certain unification schemes such as string 
theory. Whether some modified gravity theories could be nonperturbatively quantized is certainly an interesting and challenging 
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question. In this review, as an example, we first set up the s«(2)-connection dynamical formalism of STT. Then LQG method 
has been successfully extended to the STT by coupling to a polymer-like scalar field. This successful extension strongly hints 
that the nonperturbative quantization procedure might be valid even for more general modified gravity theories. At least, as 
we demonstrated, loop quantization procedure should be valid for any metric theories with a well-defined geometrical dynam- 
ics. Hence it is desirable to study the Hamiltonian formulation of modified gravity theories and try to cast those theories into 
the sM(2)-connection dynamical formalism. Then, we can naturally extend nonperturbative loop quantization method to those 
theories. 

The concrete results of this paper are summarized as follows. A general loop quantization scheme for metric modified gravity 
is first given in section [TT] Then we use STT as an example to show how our general procedure works. By doing Hamiltonian 
analysis, we have successfully derived the Hamiltonian formulation of STT of gravity from their Lagrangian formulation. The 
result shows that these theories can be naturally divided into two different sectors by the coupling parameter a>((p). In the first 
sector of u>((f>) + -3/2, the resulted canonical structure and constraint algebra of STT are similar to those of GR minimally 
coupled with a scalar field. While in the sector of co(<p) = -3/2, the feasible theories are strongly restricted and a new primary 
constraint which generating conformal transformations of spacetime is obtained. The corresponding canonical structure and 
constraint algebra are also obtained. It is worth noting that the Palatini f(H) theories are equivalent to this sector of STT. The 
successful background independent LQG relies on the key observation that GR can be cast into the s«(2)-connection dynamics. 
We have shown that the connection dynamical formalism of the STT can also be obtained by canonical transformations from the 
geometrical dynamics. Based on the connection dynamical formalism with structure group S U(2), loop quantization method 
has been successfully extended to the STT by coupling to a polymer-like scalar field. The quantum kinematical structure of 
STT is as same as that of loop quantum gravity coupled with a scalar field. Thus the important physical result that both the area 
and the volume are discrete at kinematic level remains valid for quantum STT of gravity. While the dynamics of STT is more 
general than that of LQG, the Hamiltonian constraint operators and master constraint operators for STT can also be well defined 
in both sectors respectively. In particular, in the sector «(</>) = -3/2, the extra conformal constraint can also be promoted as a 
well-defined operator. Hence the classical STT in both sectors have been successfully quantized non-perturbatively. This ensures 
the existence of the STT of gravity at fundamental quantum level. As the cosmological application of the above loop quantum 
STT, we construct a particular type of loop quantum scalar-tensor cosmology, which is the so-called Brans-Dicke cosmology. 
For simplicity, we only restrict ourselves to the sector of 10 + -|. It turns out that the classical differential equation of Brans- 
Dicke cosmology, which represents the cosmological evolution, is now replaced by quantum difference equation. The effective 
Friedmann equation of loop quantum Brans-Dicke cosmology is also given, which shows that the classical big bang singularity 
is again replaced by a quantum bounce. This effective equation lays a foundation for further phenomenological investigation to 
possible quantum gravity effects in Brans-Dicke cosmology. 

It should be noted that there are still many aspects of the connection formalism and loop quantization of modified gravity, 
which deserve discovering. Taking STT for examples, it is still desirable to derive the connection dynamics of STT by variational 
principle. The semiclassical analysis of loop quantum STT is yet to be done. In our loop quantum Brans-Dicke cosmology, some 
phenomenological issues, such as inflation, would be studied in future works. To further explore the physical contents of the 
loop quantum STT, we would also like to study its applications to black holes in future works. In addition, one would also like to 
quantize STT via the covariant spin foam approach. Furthermore, nonperturbative loop quantization of other types of modified 
gravity, such as Horava-Lifshitz theory and critical gravity etc, is also desirable. 
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